﻿ Sisteme de ordinul I şi II 1 Scopul lucrării Se studiază comportarea în domeniul timp şi frecvenţă a sistemelor de ordinul II Sisteme de ordinul I 2 Comportarea în domeniul timp a sistemelor de ordinul I Un sistem de ordinul I de tip trece-jos este descris de ecuaţia diferenţială: ()dy t (2 1) () ()yt xtτ += dt unde () irea sistemului τ se numeşte constanta de timp a xt este intrarea iar ()yt este ieş sistemului În condiţii de repaos iniţial ecuaţia (2 1) descrie un sistem liniar si invariant in timp a cărui funcţie de sistem este: 11 () , Re{}Hs s= (2 2) TJ>− τ + 1sτ Se obţin: • răspunsul la impuls: /1tτ− () ()ht e tσ = (2 3) τ • răspunsul indicial (la semnalul treaptă): t ⎡⎤ /tτ− () ( ) () (1 ) ()st h d t e tθ (2 4) θσ σ ==− ⎢⎥ ∫ 0⎣⎦ Se observă că, pe durata unei constante de timp, răspunsul la impuls scade de e ori Un sistem de ordinul I trece-sus este descris de ecuaţia diferenţială: () ()dy t dx t ()ytττ += (2 5) dt dt În condiţii de repaos iniţial, funcţia de sistem este: 1sτ (2 6) () , Re{}Hs s= TJ>− + 1sτ τ 3 Comportarea în domeniul frecvenţă a sistemelor de ordinul I Răspunsul în frecvenţă al unui sistem de ordinul I trece-jos este: 1 () ()|HHsω (2 7) == TJ TJ s jω= τ + 1jω iar al unui sistem trece-sus: jωτ == (2 8) () ()|HHsω TS TS s jω= τ 1jω+ 1Sisteme de ordinul II 4 Comportarea în domeniul timp a sistemelor de ordinul II trece-jos Ecuaţia diferenţială ce caracterizează sistemele de ordinul II trece-jos este: 2 22() ()dyt dyt ω ω++= (2 9) 2()()ytxtξω 2nnn dt dt xt este intrarea iar ()yt este ieşirea sistemului ω se numeşte pulsaţie naturala, iar ξ se unde () n numeşte gradul de amortizare al sistemului În condiţii de repaos iniţial, ecuaţia (2 9) descrie un sistem liniar si invariant în timp, a cărui funcţie de sistem este: 2 nω (2 10) ( ) Re{ }Hs sσ 022=> 2ssξω++ nn ω Polii funcţiei de sistem sunt daţi de rădăcinile ecuaţiei: 22 ++= (2 11) 20ssξω nn ω şi au expresiile: 2 −=ξωξωc (2 12) 1− nn 1+ 2 −=ξωξωc (2 13) 1− 2− nn • În cazul 01ξ , polii sunt distincţi si reali Se spune ca sistemul este în regim subcritic şi 2 1−+−=ξωξωσ nn0 În ultimele două cazuri sistemul de ordinul II este echivalent cu două sisteme de ordinul I conectate în cascadă Răspunsul la impuls al sistemului se obţine aplicând transformata Laplace inversă funcţiei (2 10) Dacă 1ξ ≠ se obţine: ct c t ⎡⎤ 12() ()ht M e e tσ =− (2 14) ⎣⎦ unde: nω (2 15) M= 2 21ξ− Pentru 1ξ = se obţine: − 2tω (2 16) () ()nht te tω = nσ Pentru 01ξ 2ssξω++ nn ω • trece-sus 2 TSAs (2 31) () Re{}Hs sσ 022TS=> 2ssξω++ nn ω • trece-bandă 2 2Asξω⋅ TB n (2 32) () Re{}Hs sσ 022TB=> 2ssξω++ nn ω • opreşte-bandă 22 ()Asω+ OB n (2 33) () Re{}Hs sσ 022OB=> 2ssξω++ nn ω Pentru valorile lui σ relaţiei (2 13) 0 vezi comentariile ce urmează 7 Desfăşurarea lucrării Vom folosi generator de semnal, • osciloscop • filtru: trece jos FTJ – led verde; trece banda FTB – led portocaliu; filtru trece sus FTS – led • rosu sursa de alimentare (±12V) pentru filtru • Iesirea 50Ω a generatorului de semnal este conectata la intrarea osciloscopului (canalul 1) si la intrarea filtrului Iesirea filtrului este vizualizata pe canalul 2 al osciloscopului 7 1 Se determină experimental caracteristicile de frecvenţă pentru trei tipuri de sisteme de ordinul doi (trece jos, trece sus, trece bandă) Se masoara caracteristica de amplitudine si de faza pentru filtre Datele iniţiale sunt: 1UV= in, amplitudinea semnalului de la intrare (sau 2V varf-la-varf) 3 4 kHzf= a de tăiere a filtrului t, frecvenţ Se completează tabelul de mai jos pentru cele trei tipuri de filtre: 6fkHz 0 1 0 4 1 4 2 4 2 9 3 2 3 4 3 6 3 9 4 4 5 4 6 4 7 9 10 20 [] [] outUV TsΔ [] Ts [] radϕ [] Defazajul se calculeaza folosind regula de trei simpla: TsΔ…… [radϕ []] Ts ……… -2π [] 2T−πΔ rad= ϕ [] T TsΔ si [Ts sunt citite folosind cursorul osciloscopului La fel si [UV (Atentie: Valorile [ ]]] out folositi VΔ 2 pentru ca semnalul de iesire este pe canalul 2!) grafic pe hârtie milimetrică caracteristicile de amplitudine si faza: 7 2 Se reprezintă [] outUV functie kf= () [] inUV rad functie kfϕ= () [] In total: 6 grafice 8 Exerciţii în Matlab Semnal sinusoidal Să se reprezinte grafic funcţiile ()sin 2 50fttπ= şi (gt f t=− Să se scrie titlul „Graficele ())() fx şi ()gx”, pe axa x se scrie „t”, iar pe axa y să se scrie ()ft şi ()gt funcţiilor () t=0: 001:0 02 f=sin(2*pi*50*t) g=-f plot(t,f,t,g,'g'),grid on title('Graficele functiilor f(t) si g(t)') xlabel('t'), ylabel('f(t) si g(t)') Exerciţiu 1⎛⎞ , pentru [ sin 2xnnπ Să se reprezinte grafic funcţia discretă: () Graficul să fie de = 0, 20n∈ ] ⎜⎟ 10⎝⎠ culoare roşie Să se scrie titlul şi identificările axelor Rprezentarea se face cu comanda stem Convoluţia semnalelor 7 12xn n n nδδ δ=+−+− şi Să se calculeze convoluţia liniară între secvenţele: [ ][][][] Se defineşte convoluţia lor liniară prin: 123hn n n n nδδ δ δ=+−+−+− [][][][][] ∞ ynxkhnk=− [] [][ ] ∑ ∞ k=− x= ; h= disp('Rezultatul convolutiei este:') y=conv(x,h) stem(y) Să se calculeze şi să se reprezinte grafic produsul de convoluţie liniară a secvenţelor: 5xn n nσσ=−−, pentru 010n≤≤ [][][] n ≤ 0, 9hn=, pentru 020n≤ () [] x=[ones(1,5),zeros(1,6)]; n=0:20; h=0 9 ^n; y=conv(x,h); subplot(2,2,1),stem(0:10,x),title('x'),grid subplot(2,2,2),stem(n,h),title('h'),grid subplot(2,1,2),stem(0:length(y)-1,y),title('y'),grid Exerciţii 1 Se dau secvenţele: 1, daca n=0,1,2,3,4,5⎧ xn= [] ⎨ 0, in rest ⎩ n+1, daca n=0,1,2⎧ hn= [] ⎨ 0, in rest ⎩ Calculaţi analitic [ yn xn hn=∗, folosind funcţia conv, şi reprezentaţi grafic rezultatul obţinut ][][] 2 Se dau secvenţele: 1, daca n=0,1,2⎧ xn= [] ⎨ 0, in rest ⎩ 5-n, daca n=0,1,2,3,4⎧ hn= [] ⎨ 0, in rest ⎩ Calculaţi analitic [ yn xn hn=∗, folosind funcţia conv, şi reprezentaţi grafic rezultatul obţinut ][][] Semnale periodice În MATLAB nu se pot genera secvenţe de lungime infinită astfel încât trebuie precizat numărul de perioade pentru o anumită secvenţă 1 [xnn= pentru 05n≤≤ (3 periode) ] 1 n=0:5; x1=n; figure(1) stem(n,x1),grid x11=[x1,x1,x1]; figure(2) 8 stem(0:(length(x11)-1),x11),grid Exerciţii Să se definească şi să se reprezinte grafic următoarele secvenţe: 1 [24xn n nσσ=−−− pentru 05n≤≤ (5 periode) ][][] 2 2 [13xn n nδδ=−−− pentru 05n≤≤ (7 periode) ][][] 3 Semnale complexe 1 Să se definească secvenţa complexă: π jn 5 pentru 20 20n−≤≤ [] 1xne= Să se reprezinte partea pară respectiv impară pentru această secvenţă n=-20:20; x1=exp(j*n*pi/5); subplot(2,1,1),stem(n,real(x1)),title('Real') subplot(2,1,2),stem(n,imag(x1)),title('Imaginar') Exerciţiu Să se definească secvenţa complexă: 2 1 [ pentru 20 20n−≤≤ − () ] 231xn n jn=− Să se reprezinte partea pară respectiv impară pentru această secvenţă Explicaţi rezultatul obţinut în urma comenzilor: plot(x1) plot(x2) 9